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We study the conditions for obtaining maximally multipartite entangled states (MMES) as non- 
degenerate eigenstates of Hamiltonians that involve only short-range interactions. We investigate 
small-size systems (with a number of qubits ranging from 3 to 5) and show some example Hamilto- 
nians with MMES as eigenstates. 
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I. INTRODUCTION 

The elusive features of multipartite entanglement are 
attracting increasing attention lately. While in the bi- 
partite case different mathematical dehnitions are phys- 
ically equivalent [lH3]j ^ unique characterization of mul- 
tipartite quantum correlations does not exist, and in- 
teresting alternative proposals are possible [IHS], that 
highlight different features of this inherently quantum 
phenomenon, including links with complexity and frus- 
tration [3]. The interest in multipartite entanglement is 
motivated by possible applications in quantum enhanced 
tasks, but also by genuine foundational aspects. 

We proposed that the multipartite entanglement of a 
system of qubits can be characterized in terms of the 
distribution function of bipartite entanglement (e.g., pu- 
rity) over all possible bipartitions of the qubits [10] . This 
led us to formulate the notion of "maximally multipar- 
tite entangled states" (MMES), as those states for which 
average purity (over all balanced bipartitions) is mini- 
mal. This notion can be extended to continuous variable 
systems [TTl [T2] and unearths novel applications, such as 
quantum teamwork |11| and controlled qubit teleporta- 
tion [13]. 

By their very definition, MMES exhibit very strong 
and distributed non-local correlations. This naturally 
leads to the following question: can MMES be obtained 
by making use of Hamiltonians that only involve local in- 
teractions? In the context of spin systems (that mostly 
concerns us here), "local" means both few-body and 
nearest-neighbors. For example, it would be remarkable 
if one could find Hamiltonians containing up to 2-body 
interaction terms, whose ground state is a MMES. If this 
were impossible, one could soften the requirement and 
ask whether can one find Hamiltonians containing up to 
2-body interaction terms, whose eigenstate is a MMES. 
This is the problem we intend to tackle in the present 
paper. A similar problem was analyzed in the context of 
GHZ states [HHH]. 

This article is organized as follows. In Section [H] we 
review the notion of maximally multipartite entangled 
state. In Section [Hi] we sketch out our general strategy 



for the search of MMES as eigenstates of Hamiltonians 
with short range interactions. In Section [TVl we explicitly 
construct Hamiltonians having a MMES as an eigenstate 
for a number of qubits ranging from 3 to 5. This is of 
interest in few-qubit applications. We conclude with a 
discussion on possible perspectives and applications. 

II. MAXIMALLY MULTIPARTITE 
ENTANGLED STATES 

Let a system of n qubits be in a pure state € "H, 
which is the only case we will consider henceforth. We 
consider a partition {A, A) of the system 5 = {1, 2, . . . n}, 
made up of and qubit, respectively, with = 
n and ua < with no loss of generality. Purity reads 

7ta = Tt{p\), (1) 

where 

PA = TrnAm{4,\) (2) 

is the reduced density matrix of party A. Purity ranges 
between 

^<^A<1, (3) 

where the upper bound 1 is reached by unentangled, fac- 
torized states with respect to the bipartition {A, A). On 
the other hand, the lower bound is obtained for maxi- 
mally bipartite entangled states, whose reduced density 
matrix is completely mixed: 

PA=^tA, (4) 

where 1a is the identity operator on the Hilbert space of 
subsystem A. 

The extension of this treatment to the multipartite sce- 
nario is based on the average purity ("potential of mul- 
tipartite entanglement" ) [TOl [T7] 
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where C"^ is the binomial coefficient, \A\ is the cardinal- 
ity of the set A and the sum is over balanced bipartitions 
TiA = ["/2], M denoting the integer part. The quantity 
ttme in Eq. (5]) measures the average bipartite entangle- 
ment over all possible balanced bipartitions and inherits 
the bounds (Sj) (with ha = [n/2]) 



2[»/2] 



<^me(|^A)) < 1- 



(6) 



A maximally multipartite entangled state (MMES) |10| 
\ip) is a minimizcr of ttme, 



TT^") =min{7rME(|^)) I m&HsAm 



(7) 



!}■ 



Given a quantum system whose (pure) state is a MMES, 
the density matrix of each one of its subsystems A d S 
is as mixed as possible (given the constraint that the to- 
tal system is in a pure state), so that the information 
contained in a MMES is as distributed as possible. The 
average purity ^ is related to the average linear en- 
tropy |17j and extends ideas put forward in O [TH] . This 
quantity has also been used to discuss generalized global 
entanglement in one-dimensional critical systems jl9ll20) . 

We can also define perfect MMES, obtained when the 
lower bound (|6| is saturated 
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We notice that a necessary and sufficient condition for a 
state to be a perfect MMES is to be maximally entangled 
with respect to balanced bipartitions. On the other hand, 
this requirement can be too strong; it can be shown that 
perfect MMES do not exist for n > 8 [T7]. 



III. GENERAL STRATEGY 



This simple observation enables us to define our prob- 
lem more precisely. The Hamiltonian ^ can be sepa- 
rated into a local part, in the sense defined above, and a 
nonlocal part, that contains all other interaction terms: 



H{e,IC) = i7ioc(e,/C) + ffnonioc(e, •'C). 



(11) 



Our desideratum is to find a set of parameters (e, /C) 
such that i?nonioc(e,'^) = 0, so that the MMES \ip) is 
a non-degenerate eigenstate (possibly the ground state). 
Clearly, this requirement might be impossible to satisfy. 
In the following we will consider some explicit examples 
for systems of 3, 4, and 5 qubits. 



IV. ENCODING A MMES INTO THE 
EIGENSTATE OF A HAMILTONIAN 

A. Three qubits 

For a system of three qubits, MMES are equivalent by 
local unitaries to the GHZ states: 
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where we have used the conventions cr^|0) — |0) and 
CT^jl) = — 11), (T^ being the third Pauli matrix. Since 
these states form a basis of the Hilbert space of the sys- 
tem, the most generic Hamiltonian can be written as 



The problem of finding a Hamiltonian involving local 
(two-body and nearest-neighbor) interactions and on-site 
external magnetic fields, one of whose eigenstates is a 
MMES is non-trivial. As a matter of fact, as we explained 
in the Introduction, MMES exhibit strongly non-local 
correlations, which, in principle, could be impossible to 
obtain by using only local terms. On the other hand, it is 
trivial to find Hamiltonians involving n-body interaction 
terms, whose ground state is an rt-qubit MMES: consider 
n qubits on a circle and the Hamiltonian 



H{e,JC) = eV + H^{IC), 



(9) 



where V = \ip){(p\ is the projection on the MMES \lp) 
and is a Hermitian operator depending on the set of 
parameters /C and satisfying 



VH^V = 0. 



(10) 



If iJ-L = and e < 0, the MMES \(p) is by construction 
the non-degenerate ground state for H. 



M, 



(16) 



where 



Vt = \Gf){Gf\ with z = 1,2, 3, 4 (17) 



are the projections on the GHZ basis states and Hm is a 
Hermitian operator containing terms of the form 

|G±)(G±| +H.C. with J. (18) 

We notice that each projection can be decomposed in two 
terms 



(19) 



where Qi contains products of two Pauli matrices, while 
Ci includes cubic terms. For example, can be decom- 



3 



posed in the following operators: 
Ci = i(|O00)(00O| + |lll)(lll|), 



'1"2 



(20) 



Ci = -(|000)(111| 



|111)(000|) 



(21) 



The decomposition in Eq. ( 19 1 enables us to rewrite the 



Hamiltonian ( 16 1 : 



H 



M- 



(22) 



We notice that cubic terms in the operators Ci are absent 
in Hm, and the C^'s are orthogonal to each other 

dCj =0 ^ j. (23) 

Thus, the only way to cancel such cubic terms is to im- 
pose 



This equality immediately implies that 
{Gt\H2\G+) = {G-\H,\Gr), 

where 



(24) 
(25) 
(26) 



denotes a local Hamiltonian (cubic couplings are absent). 
As a consequence, the state \G'^) (and any equivalent 
state by local unitaries) can never be the non-degenerate 
ground state of H2- Indeed, let us suppose that \G^) 
is the ground state of H2, whose spectrum ranges from 
Eq to -Emax- If \Gi) is also an eigenstate, then is 
degenerate, since condition (25) must hold; if \Gi) is not 
an eigenstate, then we must have 



Eo < (Gr|i?2|Gr) <E^ 



(27) 



and thus, as a consequence of Eq. (25 1, \G'{) cannot be 
the ground state. 

This result shows that it is impossible for a three- 
qubits MMES to be the non degenerate ground state of 
a local Hamiltonian. On the other hand, we now try to 
understand whether there exists a condition such that 
\G'^) is a non-degenerate (excited) eigenstate. The most 
general two-body (local) Hamiltonian is 



i<j 



By explicit calculation, it is possible to see that jG^) is 
an eigenstate of H2 if and only if the parameters satisfy 
the following conditions: 



Xij^Q for i = 1,2,3, 
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^y„=0 for i = l,2,3 

hi 
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J; 



•^13 



23 ' "2 — "^13 '^13' 
J^2, hl=K23+K32, 
Kl3, hl=Ki2+K2l. 



(29) 
(30) 

(31) 

(32) 
(33) 
(34) 



Thus, the Hamiltonian H2 has 23 free parameters. 

It is important to notice that the requirement that 
|G^) be an eigenstate does not have any influence on the 
coupling between two qubits along the z-axis; actually, 
terms of the form erf cr| are the only ones that leave the 
GHZ state invariant. 

To understand whether \Gf) is degenerate is not 
straightforward. The task can be simplified by reducing 
the number of free parameters satisfying Eqs. (29 1- (34 1. 



We consider a model in which the qubits are coupled 
along the x and z axes in a uniform external field acting 
along x: 



H 



Jk 
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(35) 



with periodic boundary conditions. The state |G^) cor- 
responds to the eigenvalue 3J, which is non-degenerate 
if and only if 



J 7^0, k^O, J^-- 



(36) 



As has been argued, the GHZ cannot be the ground state 
of Hjk'. it corresponds to the first non-degenerate excited 
state if the system is completely ferromagnetic ( J < 
and fc < 0) or if fc > and J < — fc/2. In the latter case, 
the ground state |GS) is not a MMES but has a large 
value qf_ entanglement, with ttme ^ 0.556 [Recall that 



in Eq. 



— 1/2]. Thus, we have found a range of 



values for the parameters of the local Hamiltonian such 
that the two lowest energy states contain a large amount 
of multipartite entanglement. Figure [T] shows the de- 
pendence on the parameters J and k of the potential of 
multipartite entanglement ttme for the ground state of 
the Hamiltonian (35). 



B. Four qubits 

In the case of four qubits it is known that perfect 
MMESs do not exist pBITfll^TH^ . Numerical and ana- 
lytical analyses show that the minimum of the potential 
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FIG. 1: (Color online) Potential of multipartite entanglement 
of the ground state 7rME(|GS)) for the Hamiltonian ( |35[ ). The 
plateau in the region where J and k are both positive cor- 
responds to the constant value ttme = 2/3 when the ground 
state energy is — ( J + 2k). Elsewhere, the ground state corre- 
sponds to the eigenvalue J + 2fc — 2\/ — 2Jk + 4k^. 



of multipartite entanglement is E^"^^ = 1/3 > 1/4. In 
this section we will search for local Hamiltonians having 
a non-degenerate eigenstate (possibly the ground state), 
corresponding to the uniform real state |24j 



(37) 



1=0 



determined by the coefficients 



C'^) = {1,1,1,1,1,1,-1,-1,1,-1,1,-1,-1,1,1,-1}. 

(38) 

This case can be treated in the same way as the three- 
qubit system: we can construct a basis formed by 16 
MMES, which will be labelled as |M|), all equivalent by 
local unitaries. The generic Hamiltonian acting on the 
Hilbert space of the system can be written as 
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H = Y, ^^Ml + Hm 



(39) 



1=1 



where Ai\ is the projection on the basis state |Af|). 
particular we have 



In 
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a-la^alal+ala'.a^.al). (40) 



By definition, are the expectation values on the basis 
states and Hm is a (Hermitian) linear combination of the 
mixed tensor products of the basis. Also in the case of 



four qubit, it can be shown that the MMES (37) cannot 



be a non-degenerate ground state of a local Hamiltonian. 
In fact, it can be ver ified that the three- and four-qubits 
couplings present in (40 ) and in the other projections A4\ 
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FIG. 2: (Color online) Graphic representation of all the terms 
of a Hamiltonian having IM4) as an eigenstate. The most gen- 
eral Hamiltonian is a linear combination with arbitrary cou- 
pling constants, of these terms. In this representation, solid 
lines correspond to couplings between the a;-components of 
the qubits, dotted lines between the y-components and dashed 
lines between the z-components. Dashed (red) arrows from 
qubit i to qubit j mean af cr|, dotted (blue) arrows erf aj, full 
(green) arrows cr\o^. A circle around a qubit means an in- 
teraction with an external field directed along x, a diamond 
along y, a square along z. The number "—1" inside a ring of 
4 qubits means that the two contributions in the term must 
have opposite coupling constants. The Hamiltonian (41 1 is 
made up of the the first six terms. 



are absent in Hm- As a consequence, the expectation 
value of the residual Hamiltonian on a basis state must 
be equal to that of the other three states. By an argu- 
ment similar to that used in the case of three qubits, we 
find that, if the state (37) is the ground state of a two- 



body (not necessarily local) Hamiltonian, it is at least 
four-fold degenerate. The search for the most general lo- 
cal Hamiltonian having \Ml) as an eigenstate has led to 
24 independent terms which leave the state \Ml) invari- 
ant, except for an overall constant. They are graphically 
illustrated in Fig. [2] For the meaning of the symbols, see 
the caption. 

The most evident characteristic of these terms is the 
abundance of interactions coupling Pauli matrices along 
different axis. The presence of these interactions is neces- 
sary if we want the eigenstate \Ml) to be non-degenerate. 
General conditions for non-degeneracy are very difficult 
to find if the coupling parameters are generic. Such con- 
ditions can be obtained more easily in a simplified Hamil- 
tonian, which depends only on a small number of param- 
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eters. For example, we have 



H 



Jk 
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(41) 



The Hamihonian (41) contains interactions of the form 



tjf cr| and the couphng to an external field along the z- 
axis. We notice that it does not couple qubits 3 and 4 to 
each other and, therefore, it can be implemented on an 
open chain rather than on a ring. The eigenstate \MD, 
corresponding to the eigenvalue 2 J, is non-degenerate if 
the following conditions hold 



J ^ ±V3fc. (42) 



Despite the model dependence on a small number of pa- 
rameters, it is not easy to analytically determine the po- 
sition of the excited eigenstate \M\) in the spectrum. 
Thus, we turned to a numerical analysis: after generat- 
ing 4 • 10^ random coupling parameters (uniformly dis- 
tributed in [-1,1]^), we found that the average position 
of the eigenvalue 2 J is in the centre of the energetic band; 
in the best case \M\) is the second excited state (as it 
cannot be the ground state). 



C. Five qubits 

For a system of five qubits, where perfect MMESs do 
exist [TU], we proceed in the same way as in the case 
of four qubits. We again choose to consider a uniform 
MMES [24, 



\Mk 



1 



31 



(43) 



with 



C(^) = {1,1,1,1,1,-1,-1,1,1,-1,-1,1,1,1,1,1,1,1, 
-1, -1,1, -1, 1,-1,-1, 1,-1, 1, -1,-1, 1, 1}. 

(44) 

The purity of the state jM^) is minimal for any bipar- 
tition and such a state is, therefore, a perfect MMES. 
As in the case of four qubits, we can construct a basis 
of the Hilbert space of the system which is formed only 
by MMES. The most general Hamiltonian acting on the 
Hilbert space of the system is 



32 



M, 



(45) 



where — \Ml){M^\ is a projection and Hm is an Her- 
mitian linear combination of the mixed tensor products 
of the basis vectors. 
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FIG. 3: (Color online) Graphic representation of all the el- 
ementary local Hamiltonians having \Ml) as an eigenstate. 
The meaning of the symbols is analogous to that in Fig. [2j 



In the case of five qubits the expressions are very com- 
plicated. Notwithstanding this, we can apply the same 
procedure followed in the previous subsection and draw 
some conclusions. Hamiltonians containing only two- 
body (not necessarily local) interactions have the same 
zero expectation value on each MMES of the chosen basis. 
Thus, none of them can be the non-degenerate ground 
state of such Hamiltonian. 

We found all the elementary Hamiltonians, contain- 
ing only local interactions, for which \Ml) is an eigen- 
state: they are graphically represented in Fig. [3j It 
has been found that, by combining the first nine terms 
with equal coupling constants, the eigenstate \Ml) is 
non-degenerate. Thus, we can obtain a non-degenerate 
MMES, eigenstate of a local Hamiltonian, for a system of 
five qubits without considering mixed interactions, as in 
the case of a four-qubit MMES. Moreover, we observed 
that removing any one of these nine terms is incompat- 
ible with the non-degeneracy of the MMES eigenstate. 
A numerical analysis based on a random sampling of 10^ 
sets of coupling parameters has confirmed (as in the four- 
qubits case) that the MMES jMg) is placed at the center 
of the energetic band. It turns out to be impossible to 
reach one of the low-lying excited states (we could not do 
better than placing the MMES at the 14th excited level) . 



V. CONCLUSIONS 

We have investigated whether it is possible to ob- 
tain n-qubit MMES (for 2 < n < 5) as eigenstates of 



6 



Hamiltonians involving only local (few-body and nearest- 
neighbors) interactions and fields. Since MMES exhibit 
very distributed non-local correlations, the answer to this 
problem is nontrivial. We found that MMES are the 
non-degenerate ground state only for n = 2. Already 
for n = 3 the requirement that a MMES be the non- 
degenerate ground state must be relaxed, and one finds 
that MMES can at most be the first non-degenerate ex- 
cited state. This can be interpreted as a manifestation of 
entanglement frustration [9l EH [26] . For n = 4 we found, 
in a restricted family of Hamiltonians, a MMES only as 
the second non-degenerate excited eigenstate. 

Besides its foundational interest, the present work is 
also motivated by few-qubit applications. Most, if not 
all, practically realizable quantum tasks involve only a 
very small number of qubits. The most advanced quan- 
tum applications require that these qubits be prepared 
in highly entangled states with high fidelity. One ex- 
pects that more performing applications would become 
possible by making use of MMES: some examples were 
proposed in [TTI - fT3] . This clearly calls for efficient meth- 
ods to prepare/generate MMES with large yield and ef- 
ficiency. For example, if a MMES were the ground state 
of some Hamiltonian H , then clearly one could engineer 
it by constructing H and letting the system relax toward 



its ground state. However, since this is not the case, one 
must consider a partial relaxation combined with con- 
trol techniques or devise alternative strategies. Another 
possible mechanism for generating MMES is dynamical 
rather than static. This is obviously related to the de- 
gree of complexity of a quantum circuit that generates 
MMES using, for instance, two-qubit gates. Future in- 
vestigations will be surely devoted to this subject. 

From the above mentioned point of view, due to the 
small number of qubits considered, it would be interest- 
ing to study the possibility of using presently realizable 
quantum systems with a proper engineering of the inter- 
action terms for generating MMES. Finally, the physical 
conditions and strategies that would enable one to ef- 
ficiently prepare sates with large entanglement, such as 
MMES, can be clearly generalized to other classes of en- 
tangled states. Work is in progress in this direction. 
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